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Cosmic shear as a tool for precision cosmology: minimising 
intrinsic galaxy alignment-lensing interference 
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Abstract. Cosmic shear leads to a correlation of the observed ellipticities of galaxies, an effect which is used to 
place constraints on cosmological parameters, and to explore the evolution of dark matter and dark energy in 
the universe. However, a possible systematic contaminant of the lensing signal is intrinsic galaxy alignment, with 
a correlation length of a few Mpc. Hirata & Seljak (2004) have recently demonstrated that for some models of 
intrinsic distortions, there may also be a cross-correlation between the intrinsic and lensing signals, which may 
dominate the intrinsic signal, and suppress the lensing power spectrum by several tens of percent. Unlike the 
pure intrinsic signal, this new term cannot be accounted for by neglecting or down- weighting pairs of galaxies 
which are physically close. Extending the correlation function tomography method of King & Schneider (2003) 
we illustrate how the impact of both intrinsic and cross-correlations can be significantly reduced, in the context of 
surveys with photometric redshift information. For a ground-based cosmic shear survey of ~ 100 sq. degrees with 
photometric redshifts, even in the presence of systematic contaminants at the level considered here, cosmological 
models degenerate in the f2 m — as plane can be distinguished well in excess of the 3-<r level. 
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1. Introduction 

The tidal gravitational field of mass inhomogeneities dis- 
■ torts the images of distant galaxies, resulting in corre- 
lations in their observed ellipticities. This cosmic shear 
signal depends upon cosmological parameters and the 
matter power spectrum (Blandford et al. 1991; Miralda- 
Escude 1991; Kaiser 1992). In 2000, four teams announced 
[ the first detections of cosmic shear (Bacon et al. 2000; 
Kaiser et al. 2000; van Waerbeke et al. 2000; Wittman 
et al. 2000; Maoli et al. 2001). The current and next 
generation of multi-colour surveys will cover hundreds 
or thousands of square degrees, and have the potential 
to measure the dark matter power spectrum at small 
scales (inaccessible to CMB measurements) with unprece- 
dented precision, and furthermore to probe its evolution. 
These surveys will tighten errors on cosmological param- 
eters, particularly when combined with results from the 
CMB, SNIa and galaxy surveys (e.g. van Waerbeke et al. 
2002). For example, Tereno et al. (2005) forecast the er- 
rors on parameters from CFHTLS (Canada-France-Hawaii 
Telescope Legacy Survey) coupled with CMB, primarily 
from WMAP (Wilkinson Microwave Anisotropy Probe), 
noting that many combinations benefit greatly from the 
joint constraint. 

A possible systematic contaminant of the lensing cor- 
relation function (£ L ) is intrinsic alignment, which may 



arise during the galaxy formation process. This has been 
subject to numerical, analytic and observational studies 
[e.g. Croft & Metzler 2000; Heavens et al. 2000 (HRH); 
Crittenden et al. 2001; Catelan et al. 2001; Mackey et al. 
2002; Brown et al. 2002; Jing 2002; Hui & Zhang 2002]. 
However, a great deal of uncertainty remains, with ampli- 
tude estimates spanning more than an order of magnitude. 

Recently, Hirata & Seljak (2004; hereafter HS04) 
demonstrated that there may be a finite intrinsic 
alignment-lensing cross-term (£ LI ), depending on the 
mechanism giving rise to intrinsic alignments. One such 
mechanism is tidal stretching: the shape of a galaxy is de- 
termined by the shape of a halo in which it forms - which 
is in turn related to the tidal field where the halo resides. 
If this tidal field also lenses a background galaxy, then one 
might expect a correlation between the ellipticities of the 
foreground and background galaxies. HS04 estimated that 
the cross-term can potentially suppress the lensing power 
spectrum by several tens of percent, adopting an intrinsic 
alignment model where the mean ellipticities of galaxies 
are linear functions of the tidal field in which they form. 
However, they note that this cross-term is strongly model 
dependent - for instance it vanishes for a quadratic rather 
than a linear dependence on the tidal field. 

It has been shown that photometric redshift informa- 
tion could be used to suppress by downweighting or 
ignoring galaxy pairs at approximately the same redshift 



2 



Lindsay King: Decoupling intrinsic-lensing interference 



(King & Schneider 2002), or by downweighting nearby 
pairs and subtracting a model of the intrinsic alignment 
signal from the observed ellipticity correlation function 
(Heymans & Heavens 2003). Since any cross-term does 
not arise from physically close galaxies, such schemes 
could not deal with this systematic. In King & Schneider 
(2003), hereafter KS03, a correlation function tomography 
method was developed to isolate the intrinsic and lensing- 
induced components of the galaxy ellipticity correlation 
function, assuming that photometric redshift information 
is available, so that the correlation function can be mea- 
sured between different redshift slices. 

In this paper, using the fact that the cross-term signal 
has a different dependence on redshift than either the pure 
lensing or intrinsic signals, being negligible when galax- 
ies in a pair have a small separation in redshift, we show 
that a similar procedure significantly reduces any system- 
atic suppression of the extracted lensing signal, and cor- 
responding error in the lensing analysis. 

2. General method and specific assumptions 

In this Section, the method proposed by KS03 is sum- 
marised, and extended to account for the possibility of a 
cross-term between intrinsic alignment and lensing. 



2.1. Method 

The E-mode ellipticity correlation function for galaxies 
with angular separation 9, and at true redshifts Z\, Zj, is 
composed of a lensing signal, an intrinsic signal, and a 
cross-term 

ao, Zi, zi) = £ L (0, zi, zj) + e(6, zi, Zj) + £ LI (0, z u zj) . (1) 

The origin and behaviour of £ T , and of £ LI is not yet well 
understood. However, the physics behind the pure lensing 
contribution is well established, £ L (#, Zi,Zj) being related 
to the 3-dimensional matter power spectrum P$ through 



< Z j) 



4c 4 

xR(w, Wi) R(w, Wj) 



miiiht); .IDA 





dll 



Ps 



(2tt) "\f(w) 



dw 
a 2 {w) 

1 -,w\ Jo(lO), (2) 



where Hq and f2 m are the values of the Hubble parameter 
and matter density parameter at the present epoch, and 
a(w) is the scale factor at comoving distance w, normalised 
such that a(0) = 1 today. J is the 0-th order Bessel func- 
tion of the first kind, and £ is the angular wave-vector. 
The comoving distance at Zi is denoted by Wi, and f(w) is 
the comoving angular diameter distance, which depends 
on the spatial curvature K: 



K- 1 / 2 sin(if !/2 w ) > o) 

f(w) = { w (K = 0) 

(-K)-^ 2 smh[{-K) l l 2 w] (K<0) 



(3) 



The function R(w,w') = f(w' — w)/f(w') is the ratio of 
the angular diameter distance of a source at comoving 



distance w' seen from a distance w, to that seen from 
w = 0. 

When we account for the availability of photomet- 
ric redshift estimates rather than spectroscopic ones, the 
galaxy ellipticity correlation function becomes 

£(d,z~i,Zj) = J dzi f dzjp(zi,Zj\zi,z~j,0)£(9,Zi,Zj) , (4) 

where p(zi, Zj\Ei, Zj, 9) is the probability to have true red- 
shifts Zi and Zj given photometric estimates I, and Zj. 
This is given by 

p(zi, Zj\zi,Zj,6) = (5) 

p(zi\zj) p(zj\Ej) [1 + g gg (r-)] 
/ dzi J dzj p{z l \z l )p{zj\z 3 ) [1 + £gg(r)] ' 

where £ gg is the galaxy spatial correlation function, which 
may also include redshift dependent evolution, and r is 
the comoving separation of galaxies in a pair. 

One can determine £(#, z~i, Zj) on a 3-dimcnsional grid 
of Nk angular separation bins of width A9 centred on 
9k (index K) , and Nz photometric redshift bins of width 
Az centred on each of z, (index /) and Zj (index J). We 
will refer to an observed signal using the notation £f} K , 
and to a theoretical prediction for a particular cosmol- 
ogy, intrinsic alignment mechanism, and cross-term using 

/ <rmod\ 

\? / UK' _ 

The lensing, intrinsic and cross-term correlations can 
be expressed in terms of sets of template functions A n , B n 
and C n 
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where a n , b n and c n are the amplitudes of the n-th lens- 
ing, intrinsic and cross-term template functions respec- 
tively. As noted in KS03, extra functions can be added 
as required, to span the range of plausible models. We 
describe our choice of models in Sect. 12. 31 below. A single 
index m identifies correlations between two bins with red- 
shift indices I and J, and with angular separation index 
K. There are N M = N z (N z + 1) A^/2 such independent 
measurements. The total of N = Nl + 2Vj + Nli gridded 
template models for the correlation functions can be writ- 
ten as an Nm x N design matrix M. , and their amplitudes 
as an A^-dimensional column vector Q so that 



(7) 



We recover in terms of the template functions: using 
the method of least squares, the best-fit estimates of the 
amplitudes Q n are those values Q n which minimise 



S = (e hs - MG) C- 1 (£ obs - MG) , 



(8) 
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where C is the covariance matrix. Since (£ mod ) m is a linear 
combination of template functions, the linear least squares 
estimators are 

g = (m t c^ 1 m) _1 M T c~ 1 i ohs . (9) 

With an observed correlation function £ obs , a design ma- 
trix M. containing our choice of template functions, and 
C, we can obtain Q i.e. the projection of £ obs into the 
template functions. This means that separate fits for the 
lensing, intrinsic and cross contributions can be obtained. 
Most importantly, for the purpose of using the lensing sig- 
nal in a subsequent analysis, is that we reduce the impact 
of any contamination from the intrinsic alignment-lensing 
cross-term as well as from the instrinsic alignment. 

2.2. Covariance matrix 

To evaluate 0, the covariance matrix C mm ' is required. 
Following KS03 a simplified model which neglects the cos- 
mic variance contribution will be used; thus the elements 
of C come from the intrinsic ellipticity dispersion of the 
source galaxies. This is the dominant contribution to the 
covariance at small angular scales; at larger angular scales, 
the cosmic variance terms start to dominate, with the 
transition angular scale depending on the survey geom- 
etry (Kaiser 1998; Schneider et al. 2002). 
The elements of the covariance matrix are 

2(a e 2 /2) 2 

Cram' = Tf I \ \^n,m' (1 + $1 ,I>Sj,JiSj j)\ , (10) 

where there is an extra contribution from auto-variance 
terms. Since C is diagonal, the elements of its inverse are 
simply (C~ 1 ) mm , = S mm ,/C mm . The galaxy ellipticity 
dispersion is denoted by er e . In bin m, the number of pairs 
is given by 

Af p (m)=Af F [n p(zi)Az] [n p(zj)Az] L 4 — t f - J ,(11) 

where no is the galaxy number density, p{zi) is the redshift 
probability density for redshift bin /, and L is the extent of 
the field, assumed to be square. t(9/L) is a function that 
takes into account the fact that fields have finite extent, 
and is evaluated numerically as in KS03. 

2.3. Basis models 

The simple, restricted, set of template functions used here 
is described below. 

2.3.1. Lensing basis models 

The lensing template functions A n are the gridded 
^(9, Zi, zj) for 3 models of the underlying cosmology: i) 
ACDM (f2 m = 0.3, cosmological constant density param- 
eter n A = 0.7), ii) OCDM (Q m = 0.3, 0\ = 0) and 
hi) trCDM (Q m = 1, £1a = 0); for simplicity, the power 
spectrum normalisation erg = 0.9, and shape parameter 



r = 0.21 for each model. We could choose templates close 
to ACDM, as described in Section 4. 

Starting with a 3-dimensional primordial power spec- 
trum P(k) oc k, the Bardeen et al. (1986) transfer function 
is used to describe its evolution, along with the prescrip- 
tion of Peacock & Dodds (1996) for evolution in the non- 
linear regime. The lensing correlation functions are calcu- 
lated using the relationship between the power spectrum 
and £ L given in J5J, and then integrated over the photo- 
metric redshift uncertainties as in Q). Here it is assumed 
that p(z\z) is a Gaussian with dispersion <7 p hot, centred on 
z. 

2.3.2. Intrinsic alignment basis models 

Nine template models B n for the intrinsic alignments are 
considered, the true spatial intrinsic correlation function 
being parameterised in terms of a correlation length i? C orrj 
and an exponent a: 

r){r, z) = (1 + z av ) a [exp (-r/JU)] , (12) 

where z av is the mean redshift of galaxies in a pair and r is 
their comoving separation. We use the approximation r 2 = 
(wi - w.j) 2 + 9 2 f 2 [(wi + Wj)/2\. R COII is taken to be [1, 3, 
10] hr 1 Mpc and a to be [—1, 0, 1]. The availability of pho- 
tometric redshift estimates is accounted for by integrating 
T](r, z) as in Q), and finally we obtain each of the model 
correlation functions on a grid. Note that the intrinsic 
models are calculated using the relationships for the dis- 
tances f(w) pertaining to the ACDM cosmology described 
above, and that we use £gg( r ) = (r/5 hr 1 Mpc) -18 . 

2.3.3. Cross-term basis models 

A simple toy-model prescription for the cross-term is 
adopted, parameterising £ LI directly, and evaluated for the 
ACDM cosmology. The main feature of the model is the 
dependence on angular diameter distance as in HS04; for 
fixed Zi < Zj, the cross-term at a certain angular scale 
increases with Zj. To describe the angular dependence, an 
exponential term is used so that the observed signal is 

C Zl (e,Zi,Zj) = Texp (-9/9 x ) f p(z)dz , (13) 

J Jk(w) 

where T is a normalisation and 9 X is a scale angle. The 
normalisation of the templates is fixed so that the magni- 
tude of the cross-term is 12% of the ACDM lensing signal 
in the lowest angular separation bin, with respect to red- 
shift slices at z = 0.62 and z — 1.13, and is negative 
as in HS04. Since the intrinsic alignment and cross-term 
models are currently poorly constrained, and this normal- 
isation is model dependent, it has been chosen to result 
in a cross-term that is sub-dominant to the intrinsic sig- 
nal for sources with small physical separations, and to 
dominate the intrinsic signal and have a measurable ef- 
fect on the two point correlation function for sources with 
larger physical separations - in other words a cross-term 
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that would have to be accounted for in the analysis of 
cosmic shear survey data. Note that for the templates, 
this normalisation is however arbitrary, since the signal is 
projected into the template functions, determining best- 
fit estimates of their amplitudes. Two sets of templates 
are used for the cross-term: a set of three templates with 
9 X set to [1, 5, 10] /3000 radians, and a set of nine tem- 
plates with 9 X values [0.05, 0.1, 0.3, 0.1, 1.0, 2.0, 4.0, 6.0, 
8.0] /3000 radians (where 1/3000 radian w 1.15' w 0.55 
Mpc at z = 1). Again, since predictions for the cross-term 
encompass a wide range of amplitudes, the template set 
would be expanded on as necessary. 

3. Results 

The results are presented in the context of possible future 
multi-colour cosmic shear surveys, with field size L — 14', 
hence the largest scale on which the ellipticity correlation 
function is available is v2 x 14'. A galaxy number den- 
sity of 30arcmin -2 and ellipticity dispersion of a e — 0.3 
(typical of ground-based observations) are used through- 
out. We consider surveys of various areas, the fiducial sur- 
vey having A/jr = 300 independent pointings. The value of 
Cphot = 0.1 is chosen since this is characteristic of that ob- 
tained with current SED fitting procedures such as hyper- 
z using a wide range of optical and near-infrared filters 
(Bolzonella et al. 2000). There are Nz = 65 redshift slices 
between z = 0.2 and 2.12, and Nk — 25 angular separa- 
tion bins between 0'.3 and 15'. 

The galaxy redshift distribution follows the pa- 
rameterisation from Smail et al. (1995), i.e. p(z) = 
/3/[z r /3 (3/^)](z/z ) 2 exp (-(z/z )' 3 ), where I> denotes 
the gamma function. We take (3 = 3/2 and z = 2/3 
yielding (z) sa 1. 

3.1. Minimising the contamination from the intrinsic 
and cross correlation signals 

We start with a survey oiNr — 300 fields. Our "observed" 
correlation functions comprise a lensing, an intrinsic 
and a cross contribution. The intrinsic alignment model 
for spirals from HRH, r/(r, z) = 0.012 exp(— r/1.5 /i _1 Mpc) 
was used to obtain the lensing signal £ L was calculated 
for a ACDM cosmology and £ LI for an exponential model 
with a scale angle of 4.8/3000 radians (which will be our 
fiducial cross-term model). In this subsection, the larger 
set of nine cross-term templates was used. 

Random gaussian distributed errors with dispersion 
a = C^ m were added to these correlation functions giving 
noise realisations, and best-fit parameters Q n were recov- 
ered. Fig.^and Fig.[21show the (noise-free) input and re- 
covered correlation functions between two combinations of 
redshift slices: Fig. 1 shows correlations for slices at z ~ 0.6 
and z ~ 1 and Fig. 2 shows the correlations between two 
neighbouring slices at z ~ 1. In both figures we show the 
total signal, and in Fig.Q]we also plot the lensing correla- 
tion function expected for an OCDM cosmology. 



For bins with a large separation in redshift, the in- 
trinsic signal is negligible, as expected. Similarly, for bins 
with a small separation in redshift, the cross-term is neg- 
ligible. The reduced % 2 values of the recovered fits to the 
noise realisations are « 1. Note that if we were simply to 
down-weight physically close pairs to minimise any intrin- 
sic alignment signal, that neglecting the cross-term may 
cause e.g. an underestimation of a$. 

This procedure was repeated, this time with the angu- 
lar dependence of the cross-term taken to be a power-law 
cx 9~ 15 . Note that this functional form is not contained 
in the template set, providing a check of the robustness 
of the method. Fig.|3 shows the (noise- free) input and re- 
covered correlation functions between two redshift slices 
at z ~ 0.6 and z ~ 1, along with the total signal, and the 
OCDM lensing correlation function. We also show in Fig.|3] 
the corresponding results for Af? = 3000 (i.e. a survey of 
wl60 sq. degrees, similar to the ongoing CFHTLS). Note 
than in practice, the template sets for both the IA and 
cross-term would encompass a wide range of models; this 
is particularly important for the cross-term as discussed 
below. 

3.2. Lifting the fi m — as degeneracy 

A well-known degeneracy in cosmic shear constraints is 
between tt m and erg. How accurately can a survey with 
A> = 3000, 1000 or 300 distinguish between a ACDM 
model and a near-degenerate model (fi m = 0.4, as = 0.78) 
for the lensing model, in the presence of both an intrinsic 
alignment and a cross signal? 

The gridded correlation functions for lensing in the 
ACDM and degenerate cosmologies were taken and added 
to the gridded HRH intrinsic correlation function, along 
with the fiducial cross-term model. The simulations in- 
volved generating 10000 (1000) noise realisations for the 
3000 and 1000 (300) field surveys, to obtain the_ "ob- 
served" gridded correlation functions £acdm an d £dcgcn- 
The best-fit amplitudes Q n for the template functions were 
recovered. A set of thirteen template functions was used 
containing a) the nine models for intrinsic alignments, b) 
the set of three models for the cross-term, along with c) ei- 
ther the ACDM lensing template or the degenerate model 
template in turn. We end up with four sets of best-fit am- 
plitudes for each size of survey, and corresponding values 
of \ 2 for how well the signal is modeled by projection into 
the templates: two sets generated using the ACDM (de- 
generate) lensing template, one of which is projected into 
a template set containing the ACDM lensing template, 
and the other projected into a template set containing 
the degenerate lensing template. Values were determined 
for Ax 2 (ACDM - degen) and (ii) A X 2 (degen - ACDM), 
corresponding to the difference in goodness-of-fit for the 
noise realisations when £acdm an d £dcgcn are used in the 
template set. When the fiducial ACDM model (degener- 
ate model) is the best-fit and is contained in the template 
set, values of A\ 2 should be negative. This gives a mea- 
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sure of our ability to differentiate between models using 
correlation functions between redshift slices. 

For Aff = 3000, when the model for the cosmology is 
contained in the template set, in all cases the noisy cor- 
relation functions for that cosmology are better fit. When 
Aff — 1000(300), and Iacdm i s contained in the template 
set, in 99.77 (93.9) % of cases the noisy ACDM correlation 
functions are better fit. Also, when £ degon is in the tem- 
plate set, 99.88 (94.4) % of the noisy degenerate correla- 
tion functions have better fits. Hence, within the assump- 
tions made, and in the presence of intrinsic and cross-term 
signals, these two cosmological models could be distin- 
guished well in excess of the 3-cr level when Aff = 1000, 
and just in excess of 2-cr when Aff — 300. 

For comparison, simulations were undertaken without 
any intrinsic and cross-term signal, so the only signal was 
from cosmic shear in a ACDM cosmology for a survey 
with Aff — 300. When Iacdm I s the lensing template for 
recovery, in just over 99% of cases the noisy correlation 
functions are better fit, compared with when £ dogon is the 
lensing template. The same holds for simulations where 
the cosmic shear signal is from the degenerate cosmology, 
and fdegen is the lensing template. 

Although in the foregoing example, the form of the 
cross-term is quite close to the templates, this is achiev- 
able in practice since the template set used would be large. 
Nevertheless, for the observed signal, let us next replace 
the cross-term model with exponential dependence on an- 
gular separation, with the power-law model (scale depen- 
dence oc #~ 15 ), and use the expanded set of nine exponen- 
tial cross-term templates. The ACDM model for lensing, 
HRH instrinsic alignment model, and power-law model for 
the cross-term were used to generate 10000 noise realisa- 
tions of a Aff = 1000 survey. Best- fit amplitudes for the 
template functions were determined in the same manner as 
above, first when the template set only contains a ACDM 
template, and second when it only contains the degener- 
ate template. Again, the two cosmological models can be 
distinguished in excess of the 3-cr level. 

4. Discussion and conclusions 

Recently, it was suggested by HS04 that the lensing corre- 
lation function may be contaminated, not only by intrinsic 
galaxy alignments, but also by a cross-term which results 
in a correlation between the ellipticity of a background 
lensed galaxy and a foreground galaxy. The current and 
next generation of cosmic shear surveys have the poten- 
tial to probe the evolution of the matter power spectrum, 
and to test other families of cosmological models (e.g. 
quintessence models) with similar physical manifestations 
as the concordance model (e.g. Benabed & Bernardeau 
2001; Tereno et al. 2005). In this era of precision cosmol- 
ogy, it is therefore essential to have the ability to quantify 
and isolate the contribution from both intrinsic galaxy 
alignments and the instrinsic-lensing cross-term in order 
to remove these systematics. Although the intrinsic and 
intrinsic-lensing terms are important in their own right, 




9 (arcmin) 

Fig. 1. Contributions to the two point correlation func- 
tion, as a function of angular scale, for slices at z = 0.62 
and z = 1.13. The inverted open triangles show the lensing 
correlation function for ACDM (and for comparison open 
circles show OCDM). Squares (inverted filled triangles) 
show the cross-term (intrinsic term) as described in the 
text. The sets of 10 lines show recovered lensing, cross and 
intrinsic signals for different noise realisations; the recov- 
ered instrinsic signals are indistinguishable on the figure, 
since they are very similar to each other. 
The total input signal is the curve offset below the recovered 
lensing signals. The survey has 300 fields each 14' on a side. 
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9 (arcmin) 



Fig. 2. As for Fig. 1 but for slices at z = 1.07 and z = 
1.13. The lensing correlation function for OCDM is not 
displayed. The total input signal is the uppermost curve. 



providing constraints on physical processes at the epoch 
of galaxy formation, here we have concentrated on their 
role as contaminants of the pure lensing signal. 

That the lensing, intrinsic and cross-term have differ- 
ent dependency on redshift enables the components to be 
isolated when photometric redshifts are available. In par- 
ticular, the impact of the potential cross-term suggested 
by HS04, which is impossible to remove by simply neglect- 
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(arcmin) 



Fig. 3. As for Fig. 1 but for a cross-term with a power-law 
dependence on angular separation. 




(arcmin) 

Fig. 4. As for Fig. 1 but for a cross-term with a power-law 
dependence on angular separation, and for a survey with 
3000 fields. 



ing or down-weighting physically close galaxy pairs, can 
be minimised using this type of tomographic procedure. 
Since the dependence on redshift of the lensing signal is 
more akin to the cross-term than to the intrinsic signal, 
a wide range of template functions for each of the terms, 
coupled with a large survey area, is beneficial in obtaining 
an estimate of the cross-term magnitude. 

Application to a real survey will warrant as large a set 
of template functions as possible, particularly a wide range 
of functional forms for the intrinsic and cross-terms. One 
strategy for the lensing templates could be to have cor- 
relation functions corresponding to cosmological models 
close to our best estimate for the cosmological parameters 
- perhaps finely sampling models which are degenerate in 
the f2 m — erg plane. If the reduced x 2 °f the best fit is 
significantly larger than 1, this indicates that additional 
template functions do need to be included. Additional con- 
straints on intrinsic alignment templates could come from 



direct measurements of the pure intrinsic alignment sig- 
nal in low redshift surveys where it could dominate lensing 
(e.g. Brown et al. 2002), or in surveys where spectroscopic 
redshifts are available. For the cross-term, HS04 note that 
an estimate of the density-intrinsic shear correlation has 
already been made on small scales in the context of as- 
sessing the contamination of galaxy-galaxy lensing in the 
SDSS (Hirata et al. 2004b). Such measurements on larger 
scales would help to define the template set, albeit at a 
low redshift so that the possibility of evolution would have 
to be accounted for. 

The main goal of the method developed here and in 
KS03 is to harness redshift information in order to min- 
imise the contamination of the lensing signal by any in- 
trinsic alignment or lensing-intrinsic signal. Although a 
comparison between how well discrete cosmological mod- 
els fit the data can be made (such as in the discussion 
above of models degenerate in the £l m — as plane), the 
intention is not to determine best-fit cosmological param- 
eters. In the absence of an intrinsic alignment or cross- 
term, Simon, King & Schneider (2004) have demonstrated 
how even crude redshift binning can greatly decrease er- 
rors on cosmological parameters determined by measuring 
shear correlation functions. Almost an order of magnitude 
improvement is obtained when a few redshift bins are em- 
ployed, and progressively little is gained as the number of 
bins is increased. For cosmological parameter constraint, 
a standard technique with much cruder redshift binning, 
and a much finer sampling of cosmological parameters 
could be applied to the cleaned lensing signal obtained 
from the method presented here. Such techniques are typ- 
ically analogous to the use of eq. (8), but in the form 



(14) 



where £(7r) is the theoretical shear correlation function 
for a set of cosmological parameters 7r. 

The simulations in this paper have been carried out 
for a field of 14' on a side; this is similar to the imager 
on VLT-VIMOS. Having a larger field would increase the 
number of pairs of galaxies in a particular bin from which 
the correlation function is measured, so decreasing the 
shot noise. The accuracy with which photometric redshifts 
can be determined for the galaxies in a survey also im- 
pacts on the results: upcoming multi-colour surveys such 
as the DES (Dark Energy Survey) will have redshift ac- 
curacies of dz = 0.03 — 0.07 for galaxies out to z ~ 1.1, 
depending on redshift and galaxy type (white paper at 
http://decam.fnal.gov). Cosmic shear surveys are carried 
out at optical frequencies, but when wide-field infrared 
observations of the same fields are carried out (with tele- 
scopes such as VISTA), this will greatly assist photometric 
redshift estimation. 

With our restricted set of exponential templates for the 
cross-term, and an exponential form for its true amplitude, 
a survey of ~ 60 sq. degrees, with photometric redshifts for 
the galaxies used in the analysis, would enable us to dis- 
tinguish between models degenerate in O m — erg in excess of 
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the 3-cr level. When the angular dependence of the cross- 
term is a power-law form, the same degree of distinction 
between models requires an expanded set of (exponential) 
template functions. Keeping in mind our simplified ansatz 
for the covariance matrix, this would crudely correspond 
to ~ 100 sq. degrees in practice - although given such 
a large number of fields we beat down cosmic variance. 
Accurate representations of the covariance matrix will be 
calculated using ray-tracing through large N-body simu- 
lations. 

Acknowledgements. I thank George Efstathiou, Max Pettini, 
Peter Schneider and Uros Seljak for discussions, and the ref- 
eree for a very helpful and constructive report. This work was 
supported by the Royal Society. 



van Waerbeke, L., Mellier, Y., Tereno, I. 2002, 

astro-ph/0206245 
Wittman, D.M., Tyson, J. A., Kirkman, D., Dell'Antonio, I., 

Bernstein, G. 2000, Nat, 405, 143 



References 

Bacon, D.J., Refregier, A.R., Ellis, R.S. 2000, MNRAS, 318, 
625 

Bardeen, J.M., Bond, J.R., Kaiser, N., Szalay, A.S. 1986, ApJ, 
304, 15 

Benabed, K., Bernardeau, F. 2001, Phys Rev D, 64, 8 
Blandford, R.D., Saust, A.B., Brainerd, T.G., Villumsen, J.V. 

1991, MNRAS, 251, 600 
Bolzonella, M., Miralles, J.-M., Pello, R. 2000, A&A, 363, 476 
Brown, M.L., Taylor, A.N., Hambly, N. C, Dye, S. 2002, 

MNRAS, 333, 501 
Catelan, P., Kamionkowski, M., Blandford, R.D. 2001, 

MNRAS, 320L, 7 
Crittenden, R.G., Natarajan, P., Pen, U.-L., Theuns, T. 2001, 

ApJ, 559, 552 
Croft, R., Metzler, C. 2000, ApJ, 545, 561 
Heavens, A., Refregier, A., Heymans, C. 2000, MNRAS, 319, 

649 

Heymans, C, Heavens, A. 2003, MNRAS, 339, 711 
Hirata, C, Seljak, U. 2004, Phys. Rev. D, vol. 70, Issue 6, 
063526 

Hirata, C, Mandelbaum, R., Seljak, U. et al. 2004, MNRAS, 
353, 529 

Hui, L., Zhang, J. 2002, astro-ph/0205512 
Jing, Y.P. 2002, MNRAS, 335L, 89 
Kaiser, N. 1992, ApJ, 388, 272 

Kaiser, N. 1998, ApJ, 498, 26 

Kaiser, N., Wilson, G., Luppino, G. 2000, astro-ph/0003338 
King, L., Schneider, P. 2002, A&A 396, 411 
King, L., Schneider, P. 2003, A&A 398, 23 
Mackey, J., White, M., Kamionkowski, M. 2002, MNRAS, 332, 
788 

Maoli, R., van Waerbeke, L., Mellier, Y., Schneider, P., Jain, 
B., Bernardeau, F., Erben, T., Fort, B. 2001, A&A 368, 
766 

Miralda-Escude, J. 1991, ApJ, 380, 1 

Peacock, J.A., Dodds, S.J. 1996, MNRAS, 280, L19 

Schneider, P., van Waerbeke, L., Kilbinger, M., Mellier, Y., 

2002, A&A 396, 1 
Simon, P., King, L., Schneider, P. 2004, A&A, 417, 873 
Smail, I., Hogg, D., Yan, L., Cohen, J. 1995, ApJ, 449L, 105 
Tereno,I., Dore, O., vanWaerbeke,L., Mellier,Y. 2005, A&A, 

429, 383 

van Waerbeke, L., Mellier, Y., Erben, T. et al. 2000, A&A, 358, 
30 



